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ABSTRACT 



The effects of variable properties on heat transfer in a 
film cooled turbulent boundary layer were investigated. A 
new procedure was developed to deduce adiabatic effective- 
ness from heat transfer coefficients based on the wall to 
freestream temperature difference, where both are represen- 
tative of variable property flow conditions. The new tech- 
nique was shown to be valid using data from the literature 
for injection into a turbulent boundary layer from one and 
two rows of injection holes. From these results, the varia- 
tion of the coolant to mainstream density ratio was shown to 
have a significantly greater effect on heat transfer than 
variations of viscosity and thermal conductivity. 
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NOMENCLATURE 



d 

H 

h 

h 

h 

I 

k 



f 



o 



m 

M 

P 

Pr 

q" 

Re 

T 

U 



u 

V 

W 

X 



- Injection Hole Diameter 

- Enthalpy 

- Heat Transfer Coefficient with Film Cooling 

(spanwise averaged), q'V - T^) and q"/ 

- Heat Transfer Coefficient with Film Cooling 

(spanwise averaged), q’V - T^) 

- Heat Transfer Coefficient without Film Cooling 

(spanwise averaged), q'V(T^_^ - T^) and q”/ (T^^ - T^) 

2 2 

- Injection Momentum Flux Ratio, P^-U^/p^U^ 

- Thermal Conductivity 

- Mole Fractions of Constitutant Gases 

- Blowing Ratio, p^U^/p^U^ 

- Static Pressure 

- Molecular Prandtl Number 

- Wall Heat Flux 

- Reynolds Number, 

- Absolute Temperature 

- Mean Velocity 

- Velocity Component in X Direction 

- Velocity Component in Y Direction 

- Molecular Weights of Constitutant Gases 

- Distance from Downstream Edge of Injection Holes 
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x/d - Nondimens ional Distance 

a - Thermal Diffusivity 

p - Density 

n - Adiabatic Film Cooling Effectiveness (spanwise 

averaged), eqn. 15 and 22 

y - Absolute Viscosity 

V - Kinematic Viscosity 

0 - Nondimensional Coolant Temperature, eqn. 16 

Y - Injection Angle Measured from Mainstream 

Flow Direction 

C - Nondimensional Distance Parameter, eqn. 25 

X - Second Coefficient of Viscosity 

SUBSCRIPTS 

aw - Adiabatic Wall 

c - Coolant 

e - Some Reference Value 

o - Total Condition 

r - Recovery Condition 

“ - Freestream 

SUPERSCRIPTS 

— - Averaged or Mean Value 

- Variable Property 

’ - Fluctuation from the Averaged Value 
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I. INTRODUCTION 



Over the past thirty years, gas turbines have become 
smaller and lighter with improved fuel economy. These 
desirable characteristics have been achieved by raising the 
turbine entry temperatures. To prolong useful blade life in 
this high temperature enviornment, cooling systems, thermal 
barriers, and/or different materials capable of withstanding 
higher temperatures are utilizied. 

This thesis concentrates on film cooling. Film cooling 
protects components from high gas temperatures by blanketing 
the component with a film of cooler air. This film insu- 
lates the component and acts as a heat sink. The problem 
encountered with film cooling is the prediction of the 
component's wall temperature. Two different convective heat 
transfer coefficients are currently used for this 
prediction: 



= h.CT 
f aw 



- T ) 



( 1 ) 



and 



= hr 



0«J 



- 



( 2 ) 



Here, different temperature differences are used to define h 
and h£. 

The objective of this thesis is to find a relation 
between the two heat transfer coefficients when the effects 
of variable properties are important. 
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II. BACKGROUND 



A. TURBINE BLADE COOLING FUNDMENTALS 

For a given size gas turbine, raising the turbine entry 
temperature results in a higher specific power. This 
enables the designer to operate the ’engine cycle’ at a 
higher pressure ratio, resulting in improved thermal effi- 
ciency of the thermodynamic cycle and improved overall effi- 
ciency. The net result is a lower specific fuel 
consumption. 

Raising the turbine entry temperature results in higher 
component temperatures. These components become more prone 
to failure at higher temperatures through greater oxidation 
rates of their surfaces, excessive creep, and/or by thermal 
fatigue. Figure 1 (Hennecke, 1982) shows a comparsion of 
the rising turbine entry temperature to the rising blade 
material temperature limitations. Notice that uncooled 
metal components are limited to a turbine entry temperature 
of approximately 1200 K. Current fuel economy goals require 
a turbine entry temperature of 1600 K and greater. 
Numereous cooling systems have been devised to maintain 
internal components at or below their maximum allowable 
temperature. Vanes, blades, disks, liners, and casings may 
all require some form of cooling protection. 

Figure 2 (Daniels and Schultz, 1982) illustrates two 
possible cooling medium schemes. Liquid systems use 
external sprays and/or internal convection methods. Air 
schemes use external films and/or internal convection 
methods. Liquid systems have several thermodynamic advan- 
tages over air systems. Dispite these advantages, however, 
air systems have been universally accepted due to their 
relative mechanical simplicity to implement. Three prin- 
ciple methods ‘ of air cooling currently exist. They are 
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internal convection, film, and transpiration/ effusion. 
Figure 3 (Hennecke^ 1982) illustrates these characteristic 
types of cooling for turbine blades. 

Internal convection is cooling by convection between the 
internal sections of the blade and the air. This type of 

cooling may be subdivided into several methods consisting 
of: 

(1) Straight through radial holes or channels of a single 
or multipass nature; 

(2) Enhanced heat transfer due to an enlargement of the 
internal surface area, promotion of .turbulance, or 
impingement ; 

(3) An insert that distributes the air according to the 
localized cooling requirements. 

Film cooling is localized ejection of air through 
discrete slots or holes on the blade. This results in a 
film of air around the blade which is cooler than combustion 
gas temperature, .acting as a heat sink. The film also acts 
like a thermal barrier between the hot gases and the blade's 
surface. This type of cooling may be subdivided into 
several methods consisting of: 

(1) Single films; where specific parts of the blade are 
covered by a ^Im; 

(2) Full coverage; where many films are placed close 
enough together so they cover part or all of the 
blade. 

Transpiration or effusion cooling is the bleeding of air 
through the blade at every point to form a protective film. 
The blade is made of a porous material to accomplish this 
method of cooling. 

Typically, cooling air is extracted from various stages 
in the compressor and supplied to the components that 
require cooling. Figure 4 (Hennecke, 1982) illustrates the 
coolant flow path for a single stage of a typical gas 
turbine. Typical coolant flow rates as a percentage of the 
gas turbine's total flow rate are included. Notice that the 
guide vanes and turbine blades have the highest coolant flow 
rates. ' These two components are of primary concern for they 
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have the largest amount of surface area exposed to the main- 
stream flow and consequently the highest thermal loading. 

Current applications use a combination of the internal 
convection and film cooling methods. Figure 5 (Daniels and 
Schulta, 1982) illustrates a typical guide vane cooling 
configuration. Cooling air comes in from the top and bottom 
of the vane. Notice the internal convection methods used: 

(1) Air is distributed to the internal surface; 

(2) Impingement on the surface after distribution; 

(3) Ribs and pedestals used to enlarge the surface area. 
After the air has cooled the vane by internal convection, it 
is ejected into the mainstream and -used as a film on the 
exterior surface. Initial stages utilize the internal 
convection and film cooling principles as illustrated in 
Figure 5. Latter stages may employ internal convection 
cooling solely due to the lower gas temperatures encountered 
in these stages. 

The goal in using any cooling system is to maintain the 
internal components temperatures at or below their maximum 
allowable temperature. Achievement of this goal will result 

r 

in high reliability and long service life. In addition, it 
is desirable to obtain (as closely as possible) a uniform 
blade temperature distribution to minimize thermal stresses. 

In the selection of a cooling system to be employed, two 
key limitations must be kept in mind. These limitations 
are; 



( 1 ) 

( 2 ) 



The minimal amount of cooling air should be used to 
minimize it s detrimential efrect on engine perform- 
ance ; 

The effectiveness of a particular cooling method to 
keep the component sufficiently cooled for the 
desired turbine entry temperature. 

The detrimental effect on engine performance is two 
fold. First, extracting air from the compressor reduces the 
available output power. Some work can be recovered by rein- 
troducing the cooling air back into . the mainstream flow 
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after cooling. Second, ejection of air through the blade 
has a detriment ial effect on the aerodynamic efficiency of 
the turbine. Diminished aerodynamic efficiency results from 
pressure losses and flow separation in the wakes trailing 
the blades. 

Figure 6 (Hennecke, 1982) shows the effectiveness of 
different methods used to cool components as a function of 
the turbine entry temperature. 'Cooling effectiveness' is 

defined as the ratio of actual to theoretical blade tempera- 
ture reduction. This 'cooling effectiveness’ is not only a 
function of the geometric configuration involved, but also 
the coolant flow rate. As the turbine entry temperature 
increases, so must the required cooling effectiveness. 

Figure 7 shows the cooling flow required by various 
cooling methods as a function of the turbine entry tempera- 
ture. Notice that convection cooling requires the most 
cooling flow rate while effusion requires the least at any 
given turbine entry temperature. Although this figure is 
for specific leading edge and coolant temperatures, the 
trend is the. same for other temperature" combinations. 

Figures 6 and 7 indicate that to achieve the highest 
possible turbine entry temperature with m.inimum detrimential 
effects, effusion cooling should be used. Effusion cooling, 
however, is not currently used due to structural and manu- 
facturing difficulties, and because the porous blades often 
become clogged. This is one reason why film cooling methods 
are currently employed. 

Logically, one could assume that the more coolant pumped 
into the film and the fuller the coverage, the higher the 
turbine entry temperature. However, there are limitations: 

(1) The detrimential aerodynamic effects of film cooling 
may offset the gains of a higher turbine entry 
temperature ; 

(2) "Lift-Off" may occur where the film ceases to protect 
the surface. 
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Thus, thermal design of vanes and blades with film cooling 
must be optimized. 

B. FILM COOLING 

Film cooling is the ejection of coolant into the 
boundary layer on the blade. Figure 8 shows a typical film 
cooling geometry with temperature profiles for a row of 
holes. Notice on the top view, at some distance downstream, 
the flows from adjacent jets converge to form a continuous 
film over the blade's surface. In this region, lateral 
variations become negligible, allowing a two dimensional 
analysis to be conducted with reasonable accuracy. 

The typical temperature profiles on Figure 8 illustrate 
the effect of film cooling on heat transfer. The left 
profile is without film cooling. Heat transfer is from the 
frees tream to the blade's surface. The driving force is a 
function of the freestream recovery temperature and wall 
temperature. The right profile is with film cooling. Heat 
transfer is from the freestream to the coolant and from the 
coolant to the blade's surface. The driving force is a 
function of the freestream recovery temperature, coolant 
recovery temperature, and wall temperature. As the. flow 
moves downstream from the injection site, the coolant mixes 
with the existing boundary layer. Eventually a temperature 
profile similar to the left will prevail. When this 
happens, the heat sink and insulating effects of the film 
become minimal. 

Figure 9 illustrates a typical heat transfer coefficient 
and adiabatic wall temperature distribution around a turbine 
blade. Notice that the highest temperature and heat 
transfer coefficient occurs at the nose of the blade, near 
the stagnation point, and the location of the highest 
thermal loading. Because of this, injection sites are 
located in the nose and utilize the lowest coolant tempera- 
ture. In the region around the stagnation point, the heat 
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transfer coefficients in Figure 9 are representative of a 
laminar boundary layer. At approximately one- fourth the 
distance along the suction surface, transition to a turbu- 
lent boundary layer begins, resulting in an increase in the 
adiabatic wall temperature and heat transfer coefficient 
approximately three- fourths the distance along the suction 
surface. The heat transfer coefficients along the pressure 
side in Figure 9 are representative of a turbulent boundary 
layer developing with a constant streamwise pressure 
gradient (after an initial drop from the stagnation point). 
Near the mid-chord position, the heat transfer coefficients 
begin to rise as a strong favorable pressure gradient 
develops. 

C. BOUNDARY LAYER EQUATIONS 

To predict the heat flux from the mainstream to the 
blade, the equations governing the physical situation must 
be solved. Figure 10 illustrates the coordinate system 
used. Taking a differential element and applying the 
conservation principles of mass, momentum, and energy; the 
governing equations are derived. Using the governing equa- 
tions, a solution to the heat transfer problem is consider- 
ably difficult since flow and temperature distributions are 
very complex. These equations are signif iciently simplified 
by making the following assumptions: 



(1) 


U >> V 




(2) 


3P • 

— =0 
3y 




(3) 


3T 3T 

3y 3x 




(4) 


2 

A + — U = 0 
3 




(5) 


No internal 


heat sources exist 


(6) 


No chemical 


reactions occur within the element 
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(7) Body forces are negigible 

(8) No work is performed by external fields 

Applying these assumptions reduces the governing equa- 
tions to the boundary layer equations. For turbulent 
boundary layers, the following form is for time-averaged 
velocities and variable properties: 



CONTINUITY 



MOMENTUM 



ENERGY 



where 



— (pu) + 


(pv) 


= 0 






(3) 


3x 


3y 










3"u 


9*u 


dP 


3 


3u 




(pu) — + 


(pv) — 


= . — 


+ — 


(u — - pu ' V ' ) 


(4) 


3x 


3y 


dx 


3 y 


3y 






3H 


_dP • 


3q 


dTT 




(pu) — + 


(pv) 


= u — + 




1- X 


(5) 


3x 


' 3y 


dx 


3y 


3y 




dU 












T = U 


pu ' v ' 








(6) 



3y 



and 



3T 



q = k pv’H' 

9y 



(7) 



The boundary layer equations, as presented, are highly 
nonlinear. The momentum and energy equations are coupled 
through the temperature dependency of density and viscosity. 
Thermal conductivity is also temperature dependent. The 



17 



dependence of these transport properties on temperature 
results in different velocity / temperature profiles and 
different heat transfer coefficients from those obtained by 
assuming properties which are constant through the boundary 
layer. 

Assuming constant properties exist through the boundary 
layer, the boundary layer equations may be rewritten as: 



CONTINUITY 

— -K — = 0 CS) 

3x 3y 



MOMENTUM 

3u 3"u IdP 3 3u ___ 

• u — - + V — = + (p — - u ' V ' ) (9) 

3x 3y pdx • 3y 3y 

ENERGY 

af 3T 3 3T 

V — = —(a— - v’T’) (10) 

3x 3y 3y 3y 

The momentum and energy equations ((9) and (10)) are uncou- 
pled. This means that the momentum equation can be solved 
independently of the energy equation. In addition, the 
energy equation is linear. These two items simplify the 
process of solving the boundary layer equations. 

A number of existing methods to account for variable 
properties in boundary layers without film cooling will now 
be discussed. 

Mainstream gases consist of air and the combustion 
process by-products (principally Carbon Dioxide, Carbon 
Monoxide, and water vapor). Viscosity and thermal conduc- 
tivity of this mixture can be determined though the semiem- 
pirical formula from White (1974) ; 
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n 



miMi 



mix 



i = l 



cm 



n 

T 



3 = 1 



m . 4) - - 
3 i; 



where 



Cl 



<p. . 

13 



(8 + 8W./W.)^^^ 
13 



This formula is for a mixture of n gases. W- 



C12) 

are the molec- 
The exact same 



ular weights and m^ are the mole fractions, 
formula is recommended for thermal conductivity (y is 
replaced by k). Figure 11 (White (1974)) illustrates the 
temperature variations of viscosity and thermal conductivity 
for these gases at low pressures. Notice that for a 200 R 
temperature variation, the variations in viscosity and 
thermal conductivity are not very dramatic. A power law of 
the form (Kays and Crawford (1980)): 

, 0.68 

(13) 



mix 



mix 



mix 



k. 



nix 






e e 

can be used to determine their temperature variations. For 
small temperature differences through boundary layers, say 
20 C or less, viscosity and thermal conductivity may each be 
considered to be constant. 

By considering the mainstream gases to be a perfect gas 
and from uniform static pressure through the boundary layer, 
the density variation becomes inversely proportional to the 
absolute temperature. 

^ ( 14 ) 
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Comparsion of equations 13 and 14 shows that temperature 
variations result in larger density ratio changes than u or 
k changes. 
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III. PRELIMINARY CONSIDERATIONS 



The objectives in the thermal design of gas' turbine 
blades are to obtain an acceptable stress level and design 
life. To do this, a detailed blade temperature distribution 
must be obtained. 

The thermal design of a cooled turbine blade is a 
complicated energy balance. The design process can be 
broken dovm into three major categories outlined below (in 
sequential order): 

(1) Predicting the heat flux from the hot mainstream 
gases to the blade's surface; 

(2) Performing a conduction analysis for the detailed 
mapping ox the blade s temperature distribution; 

(3) Predicting the heat flux from the blade to the 
internal cooling air. 

This thesis only addresses the first category: a relation 
is developed and verified between the two types of heat 
transfer coefficient used to predict wall temperature. 

A. CONSTANT PROPERTY CONSIDERATIONS 

Equations 1 and 2 are used to calculate the blade's wall 
temperature after the heat flux is obtained by solving the 
boundary layer equations. 



q" = h.(T - T ) 
^ aw 

and 

q” = h(T - T ) 



Cl) 

C) 



Both equations reflect the presence of film cooling. 

In developing thermal boundary layers with film cooling, 
an adiabatic wall boundary condition can be used to describe 
the convection process. Using this concept reduces the 
driving potential from a three temperature potential problem 
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to a two temperature potential problem. Equation 1 repre- 
sents this concept. Here, h£ is an isoenergtic heat 
transfer coefficient. The adiabatic wall temperature, T , 
is the wall temperature when there is no heat transfer to 
the wall. Adiabatic wall temperature is usually expressed 
nondimens icmally in terms of the adiabatic film cooling 
effectiveness, nJ 



n 



aw 



roo 



T* 'T 

1 • i. 

rc r® 



CIS) 



Since the entire surface is considered adiabatic, T is a 

aw 

global parameter. As a result, equation 1 is dependent on 
the entire convection process and the thermal boundary 
layer. 

Compared to the large temperature differences between 
the mainstream gases and the blade's surface, the surface 
temperature distribution can be considered relatively isoth- 
ermal. An isothermal wall is ideal since thermal stresses 
are minimal. Equation 2 represents this concept. Here, h 
is defined by the freestream stagnation temperature and the 
wall temperature difference for a speciJEic coolant tempera- 
ture. The coolant temperature is usually expressed nondi- 
mensionally as 9: 

T - T 
rc ■ r® 

9 = (16) 



There has been considerable controversy over these two 
methods of predicting the wall's temperature. Two key 
issues of this controversy deal with how the temperature 
difference is defined and the effects of changing coolant 
parameters (temperature, density, mass flow rate, etc.). In 
equation 1, when the temperature difference equals zero, 
there is no heat transfer. This is not necessarily true if 
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equation 2 is used. To account for viscous dissipation, the 
freestream stagnation temperature, Tq„ , is replaced by the 
freestream recovery temperature, heat transfer 
coefficient, h, is then redefined as: 



= hC 



r® 



• T ) 
w 



n") 



Changing the film cooling parameters will have different 
effects on these equations. Generally, the recovery temper- 
ature is determined from freestream conditions (temperature, 
velocity, density), however, if it is based on near wall 
conditions , the recovery temperature in a film cooled 
boundary layer is the same as the adiabatic wall tempera- 
ture. In equation 1, both h^ and are affected, whereas, 
in equation 17 only h is affected. 

For a given location, the heat flux will be the same for 
either method. As such, equation 1 is set equal to equation 
17 and a relation between the heat transfer coefficients 
exists. Using the definitions of ri and 9, and dividing by 
the heat transfer coefficient for the case of no film 
cooling, then yields: 




This equation relates adiabatic and isothermal wall boundary 
conditions for properties that are constant through the film 
cooled boundary layer. This relationship is plotted on 
Figure 12 (Ville, Cunat , and Prichards, 1978). Note the 
linearity between h/h^ and 0. This linearity is due to the 
fact that the energy equation is linear for constant proper- 
ties. This equation is valid for a given dow’nstream loca- 
tion, w’here adiabatic conditions are represented by the 
horisontal intercept, and isothermal conditions by the line. 
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On Figure 12-, the axis intercept points have special 
meaning. The h/hQ axis intercept point represents the heat 
transfer coefficient when the coolant is injected at the 
same temperature as the mainstream (T^j-=T^oo» 0=0). The 9 
axis intercept point indicates a wall condition that is 
locally adiabatic. At this point, 9 equals the inverse of 
the effectiveness, 1. 

Figure 12 is only good for one downstream location. As 
you move further downstream, the slope of the curve becomes 
less and may eventually become horziontal indicating no 
cooling effect from the film. This is a result of a rising 
film temperature due to heat transfer, diffusion, and mixing 
with mainstream gases. As the 'film temperature rises, the 
value of 0 decreases and the value of n increases. 



B. VARIABLE PROPERTY CONSIDERATIONS 

With variable properties, the momentum and energy 
boundary layer equations become coupled because of the 
density term. In addition, the energy equation becomes 
highly nonlinear. Because of this, the linearity exhibited 
on Figure 12 no longer exists. Ligrani and Camci (1985) 
illustrated this nonlinearity that exists between h/h^ and 
9. They also derived an empirical relation to relate the 
constant property data of Jabbari and Goldstein (1978) to 
their variable property data. 

To find the relation between h/h^^ and 9 for the variable 
property case, we start with the defining equations for the 
heat transfer coefficients. 

For variable properties through the boundary layer, the 
variable property isoenergetic heat transfer coeffi- 
cient, h^, is defined using: 

■ hf'Lw - TJ (19) 



The heat transfer coefficient, h, is then defined using: 



q” = h(T 
n T' 



T ) 



( 20 ) 
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The heat flux is now obtained by solving the variable prop- 
erty boundary layer equations. At a given location, the 
heat flux from equations 19 and 20 will be the same. 
Equating these gives: 

- T,,) = h(T - T ) C-1) 

r aw 



Now the adiabatic wall temperature is expressed in terms of 
an adiabatic film cooling effectiveness defined as: 



n 




T . T 

‘re 'r® 



( 22 ) 



9 is defined as before (equation 16). By substitution and 
rearrangement the variable property relation between h/h^ 
and 9 then becomes: 



h 



h. 





(23) 



This equation is for one downstream location at one 
value of the density ratio, P^/P^^, and one value of any of 
the following blowing ratio, M; momentum flux ratio, I; or 
the velocity ratio, u_/u„. A different line results from 
equation 23 for each value of the effectiveness, n , and the 
normalized heat transfer coefficient, h£/hQ. The effective- 
ness, ri, may have a strong dependence on the density ratio, 
P^/ Pgg, as indicated by the single row film cooling data of 
Pedersen, Eckect , and Goldstein (1977). 

Figure 13 is a domain map showing density ratio verses 
blowing ratio of the Pedersen, Eckect, and Goldstein (1977) 
data. The box symbols indicate that effectiveness has 
little or no dependence on the density ratio for low blowing 
ratios. The circle symbols indicate that effectiveness has 
a strong dependence on the density ratio. The triangles 
indicate that the effect of the density ratio on the effec- 
tiveness and heat transfer is highly nonlinear. 
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Figure 14 shows one possible effect of variable proper- 
ties on h/h^ verses 0 data. Here, each of the three linear 
curves represents equation 23 with a different n. The same 
h£/ho was. used for each line to simplify Figure 14. Now, 
each line is also for one value of the density ratio, /p„ , 
which is equivalent to T^^/Tc. The variation of the density 
ratio with 9 is then given by: 



P / P 



T /T 
w r® 




(24) 



if T../T_ = T /T for constant T . The 9 for a given 
density ratio will then be a straight vertical line 
eminating from the horizontal axis. The intersection of 
this vertical line and equation 23 will be where a given 
h/hj^ data point is located. The curved line represents the 
collection of the h/l^ verses 9 data points. It curves 
downward because of the way the data collects. As the 
density ratio increases (vertical lines move to the right), 
the effectiveness increases.- This causes the 0 axis inter- 
cept point (1/q) to move to the left. 
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IV. RESULTS 



Using equation 23, the adiabatic effectiveness was used 
to deduce using the data of Pedersen, Eckect, and 
Goldstein (1977) for one row of injection holes. In addi- 

/N 

tion, the adiabatic effectiveness was deduced from h/h and 
h£/h^ from the data of Ligrani and Camci (1985) for two rows 
of injection holes. 

In the following sections, the results are presented. 

A. DATA OF PEDERSEN, ECKECT, AND GOLDSTEIN (1977) 

Figure 15 shows the Pedersen, Eckect, and Goldstein 

/V 

(1977) data used to deduce h/h^. Figure 13 gives the 
blowing ratio and the density ratio values for which n was 
selected for deducing h/h^. 

A value of 1. 0 was assumed for the h/h^ axis intercept. 
Although in many cases, a reasonable assumption, Metzger, 
Carper, and Swank (1968) showed that this axis intercept may 
not be equal to 1.0. They showed that the injection has an 
effect on the heat transfer coefficient even when the 
coolant is injected at the frees tream temperature, primarily 
due to boundary layer thicking. 

Normalized heat transfer coefficients as a function of 9 
are shown on Figures 16 thru 19 for the Pedersen, Eckect, 
and Goldstein (1977) data. Four axes are used to denote 
four downstream locations (x/d). The light linear lines 
represent equation 23 for a different value of effective- 
ness, ri- The dark lines represent the collection of h/h,^ 
verses 0 data points. 

Figure 16 is for a blowing ratio, M, of 0.213. This 
figure shows that h/h^ verses 0 is steeper at small x/d and 
becomes flatter as the distance from the injection site 
increases (increasing x/d). This figure shows reasonable 
linearity with some nonlinearity near the injection site. 
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Figure 17 is for a blowing ratio, M, of 0.515. This 
figure shows the same trends as with Figure 16. Only now, 
the nonlinearity near the injection site is more pronounced 
and is also evident at a location further downstream. 

Figure 18 is for a blowing ratio, M, of 1.05. This 

figure shows that the variation of h/h^ with 9 is nonlinear. 

As the distance from the injection site increases 
(increasing x/d), the h/h^ variation with 9 becomes more 
horizontal. Notice that, in general, at a particular 
blovring ratio, the slope becomes more negative as 9 

increases. 

Figure 19 is for a blowing ratio, M, of 1.96. This 

figure shows a highly nonlinear behavior due to the ’lift- 
off' of film cooling jets from the wall, and possible pene- 
tration through the boundary layer. 

B. DATA OF LIGRANI AND CAMCI (1985) 

Figure 20 ( Ligrani and Camci, 1985) illustrates the 

downward sloping trend for variable property h/h^ verses 9 
data. This figure also shows the linear h/h^ verses 9 data 
(constant property case). Notice the similarity between the 
variable property trend here and the trends of the Pedersen, 
Eckect, and Goldstein (1977) data on Figure 18. 

^ /N 

Using h/h^ and h^/h^ from Ligrani and Camci ( 1985), n 
was deduced from equation 23. Results are plotted as a 
function of x/d in Figures 21 thru 31, and as a function 
of C in Figures 32 thru 42. The later function, is a 

downstream distance normalized by the blowing ratio. It is 
defined by: 



C 



X + 1 . 909<1 u n ^ 

■ (i^c 

MS u 



( 25 ) 



where S = 7id/6 is the width of an equivalent two dimensional 
slot (Jabbari and Goldstein, 1978) . Goldstein (1971) gives 
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two equations to -predict laterally averaged effective- 
ness, n> for two dimensional film cooling. These 
predictions are given by: 



fT = (1 + 0.2490 



- 0.8 



and 



1.9Pr 



2/3 



b ~ DR 

1 + 0.329C B 



where 



6 



1 ^ l.SdO*"^) (Re^(-£))sinT 

u_ 



(26) 



(27) 



(28) 



These predictions are plotted on the n verses ? graphs 
(Figures 32 thru 42). 

Figures 21 thru 42 are for different blowing ratios (M = 
1.242, 1.003, 0.481, 0.661, 1.557), momentum flux ratios (I 
= 1.382, 0.925, 0.417), and velocity ratios (u^/u^ = 0.630, 
1.113, 0.922). In general, these figures show that as the 
blowing ratio, momentum flux ratio, and velocity ratio 
increase, the adiabatic film cooling effectiveness 
increases. The largest increase in n occurs for the highest 
value of 0. This corresponds to the case of the highest 
density ration. On. Figures 21 thru 31, the crossover point 
between the 01 and 02 curves moves downstream (from x/d = 20 
to x/d = 40) as the momentum flux ratio and velocity ratio 
increase. On Figures 32 thru 42, as the blowing ratio, 
momentum flux ratio, and the velocity ratio decrease, the 
data moves to a higher value and spreads out. All data 
points on the n verses K graphs are below the predictions of 
equations 26 and 27 but do follow the same trend as the 
predictions. 
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V. SUMMARY 



A relation was derived that relates the normalized heat 
transfer coefficients h and h£. These heat transfer coeffi- 
cients are for isothermal wall boundary conditions, h, and 
adiabatic wall boundary conditions. h£. Each are for when 
variable properties are considered to exist through the 
boundary layer. This relation was presented as: 



AH Ah 




(2.-) 



When the effects of variable properties are important in 
film cooled turbulent boundary layers, the variable property 

/N 

effectiveness, n , and the above equation are valid for one 
density ratio and blowing ratio only. 

The normalized heat transfer coefficient, h/hQ , was 
deduced from the Pedersen, Eckect, and Goldstein (1977) data 
on injection with one row of holes. The adiabatic film 

/N 

cooling effectiveness, n, was deduced from the Ligrani and 
Camci (1985) data on injection with two rows of holes. The 
data deduced from equation 23 showed the nonlinear trend of 

/N 

h/h_ verses 9 and as the density ratio, p_/p_, varied. 

The effect of density ratio on h/ho verses 0 results in 
a change of the characteristic curve from linear to nonli- 
near. This change is not always evident. It's existance or 
extent is dependent upon the injection parameters (blowing 
ratio, momentum flux ratio, and velocity ratio) and the 
downstream distance. This effect is due to the nonlinearity 
introduced into the energy equation and the coupliing of the 
momentum and energy equations when variable properties 
through the boundary layer are considered. 
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VI. RECOMMENDATIONS 



Experiments need to be conducted to determine the accu- 
racy and dependability of the proposed variable property 
relation on injection system geometry, injection flow rate 
parameters, fluid mechanics of the flow, the thermal 
boundary layer conditions, and other parameters. In order 
to completely verify the proposed equation, simultaneous 
measurements of variable property effectiveness and variable 
property heat transfer coefficients (based on the wall and 
recovery frees tream temperature difference) need to be 
conducted. 
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APPENDIX 

FIGURES 
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Fig. 1 The Development of the Turbine Entry Temperature. 
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Fig. 2 Methods of Turbine Blade Cooling. 
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Fig. 3 Characteristic Types of Blade 
Cooling Configurations. 
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Fig. 4 Coolant Flow in a Typical 
High Pressure Turbine Stage. 
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Fig. 5 Coolant Flow in a Typical 
Nozzle Guide Vane. 
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Fig. 6 Turbine Entry Temperature Verses 
Required Cooling Effectiveness, 
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Fig. 7 Comparative Cooling Requirements. 
Leading Edge Temperature 1000 C. 
Cooling Air Temperature 260 C. 
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Fig. 8 Film Cooling by a Row of Holes 
with Sample Temperature Profiles. 
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Fig. 9 Distribution of Heat Transfer Coefficient and Adaibatic Wall 

Temperature Around a Typical Blade. 




Fig. 10 Curvilinear System of Coordinates. 
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Fig. 11 Viscosity and Thermal Conductivity 
of Common Gases at Low Pressures. 
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Fig. 12 Linear Variation of Heat Transfer 
Coefficient with Nondimens ional 
Coolant Temperature. 
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Fig. 13 Domain Map of Blowing Ratio and Density 
Ratio for Pedersen, Eckect, and Goldstein (1977) Data. 
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Fig. 14 Nonlinear Variation of h/h with 0. 
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Fig. 15 Effectiveness Verses x/d from 
Pedersen, Eckect, and Goldstein (1977). 
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Variation of li/h with 0 and x/d for 
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Fig. 17 Variation of H/Iiq with 0 and x/d for M-0.515. 
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Fig. 18 Variation of Ii/Iiq with 0 and x/d for 
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Fig. 19 Variatioii of witli 0 and x/d for 




Fig. 20 Variation of h/h with 0 for Two Rows 
of Injection from Ligrani and Camci (1985). 
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Fig. 21 Variation of n witli x/d for M=1.2A2, 
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Fig. 22 Variation of fj witli x/cl for M=1.003. 



o 




55 



Fig. 23 Variation of 
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Fig. 24 Variation of 
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Fig. 25 Variation of n with x/cl for H=1.557. 
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Fig. 26 Variation of n witli x/d for 1=1. 382. 
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Fig. 27 Variation of n with x/d for 1=0.925. 
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Variation of n with x/d for 1=0.417. 
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Fig. 29 Variation of n with x/d for u /u =0.630. 
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Fig. 30 Variation of n witli x/d for /u^ = 1. 113. 
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Fig. 31 Variation of p with x/d for /u^„ = 0. 922. 
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Fig. 32 Variation of p with f for M=1.242. 
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Fig. 33 Variation of n with ^ for M=1.003. 
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Fig. 34 Variation of n with 
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Fig. 35 Variation of n with f| for H=0.661. 
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Fig. 36 Variation of n with C for M=1.557. 
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Fig. 37 Variation of ^ with ? for 1=1.382. 
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Fig. 38 Variation of n with C for 1=0.925. 



o 




X 



< 



71 



Fig. 39 Variation of n with C for 1=0. A17. 
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Fig. 40 Variation of n with ^ for u /u =0. 630, 
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Fig. Variation of n with C for u /u =1. 113. 
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Fig. 42 Variation of ri with ^ for u^/u =0.922. 
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